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(57) Abstract 

Method and apparatus for 
determining the instantaneous 
location, the orientation of 
an object moving through a 
three-dimensional space by 
applying to the object a coil 
assembly including a plurality of 
sensor coils (20) having axes of 
known orientation with respect to 
each other including components 
in the three orthogonal planes; 
generating a time-varying, 
three-dimensional magnetic 
field gradient having known 
instantaneous values of magnitude 
and direction; applying the 
magnetic field gradient to the space, 
and object moving therethrough 
to induce electrical potentials in 
the sensor coils; measuring the 
instantaneous values of the induced 
electrical potentials generated in 
the sensor coils; processing the 
measured instantaneous values 
generated in the sensor coils 
together with the known magnitude, 
direction of the generated magnetic 
field gradient, the known relative 
orientation of the sensor coils in 
the coil assembly to compute the instantaneous 
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METHOD AND APPARATUS TO ESTIMATE LOCATION AND ORIENTATION OF 
OBJECTS DURING MAGNETIC RESONANCE IMAGING 

Field of the Invention 

5 

The invention relates to methodology and apparatus to determine the location and orientation of an 
object, for example a medical device, located inside or outside a body, while the body is being 
scanned by magnetic resonance imaging (MRI). More specifically, the invention enables estimation 
of the location and orientation of various devices(e.g. catheters, surgery instruments, biopsy 
1 0 needles, etc.) by measuring voltages induced by time-variable magnetic fields in a set of miniature 
coils. Such time-variable magnetic fields are generated by an MRI scanner during its normal 
operation. 

Background of the Invention 

15 

Minimally invasive procedures: Minimally-invasive diagnostic or interventional procedures require 
either direct visual viewing or indirect imaging of the field of operation and determination of the 
location and orientation of the operational device. For example, laparoscopic interventions are 
controlled by direct viewing of the operational field with rigid endoscopes, while flexible 

20 endoscopes are commonly used for diagnostic and interventional procedures within the gastro- 
intestinal tract. Vascular catheters are manipulated and manoeuvred by the operator, with real-time 
X-ray imaging to present the catheter location and orientation. Ultrasound imaging and new real- 
time MRI and CT scanners are used to guide diagnostic procedures (e.g. aspiration and biopsy) and 
therapeutic interventions (e.g. ablation, local drug delivery) with deep targets. While the previous 

25 examples provide either direct (optical) or indirect (imaging) view of the operation field and the 
device, another approach is based on remote sensing of the device with mechanical, optical or 
electromagnetic means to determine the location and orientation of the device inside the body. 

Stereotaxis: Computer-assisted stereotaxis is a valuable technique for performing diagnostic and 
30 interventional procedures, most typically with the brain. The concept behind the technique is to 
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have real-time measurement of the device location in the same coordinate system as an image of 
the field of operation. The current location of the device and its future path are presented in real- 
time on the image and provide the operator with feed-back to manipulate the device with minimal 
damage to the organs. During traditional stereotaxis, the patient wears a special halo-like 
5 headframe, which provides the common coordinate system, and CT or MRI scans are performed 
to create a three-dimensional computer image that provides the exact location of the target (e.g. 
tumour) in relation to the headframe. The device is mechanically attached to the frame and sensors 
provide its location in relation to the head frame. When this technique is used for biopsy or 
minimally-invasive surgery of the brain, it guides the surgeon in determining where to make a small 

10 hole in the skull to reach the target Newer technology is the frameless technique, using a 
navigational wand without the headframe (e.g. Nitin Patel and David Sandeman, "A Simple 
Trajectory Guidance Device that Assists Freehand and Interactive Image Guided Biopsy of Small 
Deep Intracranial Targets", Comp Aid Surg 2:186-192, 1997). In this technique remote sensing 
system (e.g. light sources and sensors) provides the real-time location of the device with respect to 

15 the image coordinate system. Yet both the stereotaxis and the frameless techniques are typically 
limited to the use of rigid devices like needles or biopsy forceps since their adequate operation 
requires either mechanical attachments or line of sight between the light sources and the sensors. 

Electrom agnetic remote sensing: Newer remote sensing techniques are based on electromagnetism. 

20 For example, Acker et al (US Patent 5,558,091) disclose such a method and apparatus to determine 
the position and orientation of a device inside the body. This method uses magnetic fields generated 
by Helmholtz coils, and a set of orthogonal sensors to measure components of these fields and to 
determine the position and orientation from these measurements. The measurement of the magnetic 
field components is based on Hall effect and requires exciting currents in the sensors in order to 

25 generate the measured signals. The technique requires control of the external magnetic fields and 
either steady-state or oscillating fields, for the induced voltages to reach a state of equilibrium. 
These requirements prevent, or greatly complicate, the use of this technique with magnetic fields 
generated by the MRI system, and requires the addition of a dedicated set of coils to generate the 
required magnetic fields. 



30 
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A different approach for remote sensing of location is disclosed by Pfeiler et al. (US Patent 
5,042,486) and is further used by Ben-Haim for intra-body mapping (US Patent 5,391,199). Their 
technology is based on generating weak radio-frequency (RF) signals from three different 
transmitters, receiving the signals through an RF antenna inside the device, and calculating the 
5 distances from the transmitters, which define the spatial location of the device. As with the previous 
methodology, the application of the technology to MRI is problematic due to the simultaneous use 
of RF signals by the MR scanning. Potential difficulties are the heating of the receiving antenna in 
the device by the high amplitude excitation RF transmissions of the MRI scanner and artifacts in 
the MR image. 

10 

Dumoulin and colleagues disclose another approach to determine the location of a device, 
using a small receiving coil which is sensitive to near-neighbourhood emitted RF signal during the 
MR imaging process (Dumoulin CL, Darro RD, Souza SP, "Magnetic resonance tracking", in 
Interventional MR, edited by Jolesz FA and Young IY, Mosby, 1 998). This method cannot directly 
1 5 determine the orientation of the device, and may be subject to similar difficulties as the previous 
technology, including heating of the coil. 

Interventional MRI: Many of the advantages of MRI that make it a powerful clinical imaging tool 
are also valuable during interventional procedures. The lack of ionizing radiation and the oblique 

20 and multi-planar imaging capabilities are particularly useful during invasive procedures. The 
absence of beam-hardening artifacts from bone allows complex approaches to anatomic regions that 
may be difficult or impossible with other imaging techniques such as conventional CT. Perhaps the 
greatest advantage of MRI is the superior soft-tissue contrast resolution, which allows early and 
sensitive detection of tissue changes during interventional procedures. Many experts now consider 

25 MRI to be one of the most powerful imaging techniques to guide interventional interstitial 
procedures, and in some cases even endo vascular or endoluminal procedures (Y oshimi Anzai, Rex 
Hamilton, Shantanu Sinha, Antonio DeSalles, Keith Black, Robert Lufkin, "Interventional MRI for 
Head and Neck Cancer and Other Applications", Advances in Oncology, May 1995, Vol 1 1 No. 
2). 

30 
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From the presented background on current methodologies, one can define the ideal system 
for minimal invasive procedures: It should provide real-time, 3-dimensional, non-ionizing imaging 
(like MRI or ultrasound) as feed-back to the user for optimal insertion and intervention; it should 
implement flexible, miniaturized devices which are remotely sensed to provide their location and 
5 orientation. By combining a composite image of the field of operation and the device location and 
orientation, the operator can navigate and manipulate the device without direct vision of the field 
of operation and the device. This may facilitate the use of minimal invasive intervention in the brain 
or other organs. 

10 Objects and Summary of the Invention 

An object of the present invention is to provide a novel method and apparatus for 
determining the instantaneous location and orientation of an object moving through a three- 
dimensional space, which method and apparatus have advantages in one or more of the above 
15 respects. 

Another object of the present invention is to provide such a method and apparatus which 
is particularly useful in MRI systems by making use of a basic universal component of the MRI 
system, namely the time-varying magnetic gradients which are typically generated by a set of three 
20 orthogonal electromagnetic coils in such systems. 



determining the instantaneous location and orientation of an object moving through a three- 
dimensional space, comprising: applying to the object a coil assembly including a plurality of 

25 sensor coils having axes of known orientation with respect to each other and including components 
in the three orthogonal planes; generating a time-varying, three-dimensional magnetic field 
gradient having known instantaneous values of magnitude and direction; applying the magnetic 
field gradient to the space and the object moving therethrough to induce electrical potentials in the 
sensor coils; measuring the instantaneous values of the induced electrical potentials generated in 

30 the sensor coils; and processing the measured instantaneous values generated in the sensor coils, 
together with the known magnitude and direction of the generated magnetic field gradient and the 



According to one aspect of the present invention, there is provided a method of 
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known relative orientation of the sensor coils in the coil assembly, to compute the instantaneous 
location and orientation of the object within the space. 

The above-described method is particularly useful in MRI systems, wherein the magnetic 
5 field gradient is generated by activating the gradient coils of an MRI scanner, and the invention is 
therefore described below with respect to such a system. 

According to further features in the described preferred embodiment, the magnetic field 
gradient is generated by activating three orthogonally disposed pairs of gradient coils according to 
10 a predetermined activating pattern; and the measured instantaneous values of the induced electrical 
potentials generated in the sensor coils are processed, together with the predetermined activating 
pattern of the gradient coils and the known relative orientation of the sensor coils, to provide an 
estimate of the location and orientation of the object. 

15 According to another aspect of the present invention, there is provided apparatus for 

determining the instantaneous location and orientation of an object moving through a three- 
dimensional space, comprising: a coil assembly carried by the object and including a plurality of 
sensor coils having axes of known orientation with respect to each other and including components 
in the three orthogonal planes; a magnetic field generator for generating a time-varying, three- 

20 dimensional magnetic field gradient having known instantaneous values of magnitude and direction 
in the space and the object moving therethrough to induce electrical potentials in the sensor coils; 
means for measuring the instantaneous values of the induced electrical potentials generated in the 
sensor coils; and a processor for processing the measured instantaneous values generated in the 
sensor coils, together with the known magnitude and direction of the generated magnetic field 

25 gradient and the known relative orientation of the sensor coils in the coil assembly, to compute the 
instantaneous location and orientation of said object within said space. 

The disclosed methodology and apparatus enable the estimation of the location and 
orientation of an object or a device by using a set of miniature, preferably (but not necessarily) 
30 orthogonal coils. The simplest, preferred embodiment has a set of three orthogonal coils. However 
more complex coil sets, for example a set of three orthogonal pairs of parallel coils, can improve 
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the accuracy of the tracking with a higher cost of the system. To simplify the presentation, the 
following disclosure deals specifically with a set of three orthogonal coils, and also refers to the 
more complex configuration of three orthogonal pairs of coils. However the same concepts can be 
applied to various combinations of coils by anyone familiar with the field of the invention. 

Hie time change of magnetic flux through a coil induces electromotive force (i.e. electric 
potential) across the coil (Faraday Law of electromagnetism). MRI scanners generate time-variable 
magnetic fields to create magnetic gradients in the scanned volume. By measuring the induced 
electric potentials in the three orthogonal coils (or pairs of coils), and by getting the time pattern 
10 of the generated magnetic gradients as input from the MRI scanner, both the location and 
orientation of the device can be estimated. 

The present invention has significant advantages over existing methodologies. Compared 
with stereotaxis, either the frame or frameless techniques, the new methodology enables the use of 

1 5 devices like catheters or surgical instrumentation without the need for direct line of sight with the 
device. Unlike the remote electromagnetic localization methodology of Acker et al the present 
invention is based on measurement of voltages induced by a set of time-varying electromagnetic 
gradient fields in a set of coils (Faraday Law), rather than the need to use homogenous and gradient 
fields which induce voltages in a set of miniature conductors carrying electrical current (Hall 

20 effect). Thus, the present invention is totally passive, it does not require any excitation of the 
sensors, nor the use of dedicated magnetic fields, and the requirement for time-variable magnetic 
fields is satisfied with virtually any MRI scanning protocol which is in routine clinical use. The 
methods disclosed by Pfeiler et al and Dumouline et al require the use of two sensors to measure 
orientations and thus have limited accuracy of orientation estimation, while the present invention 

25 uses a sensor which provides simultaneously accurate orientation and location tracking. Unlike 
existing optical tracking systems, there is no limitation on the number of sensors being used, and 
there is no need to maintain a line of sight between the sensor and the tracking apparatus. All other 
tracking methodologies are based on their own reference system, and should be aligned with the 
MRI coordinate system by a time-consuming registration procedure. The disclosed tracking 

30 methodology does not require registration since it uses the same set of gradient coils which are used 
by the MRI scanner for spatial encoding of the images. 



• 
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Brief Description of the Drawings 

Figure 1 A provides a block diagram description of the invented apparatus which includes 
a processing and control unit (10), a sensor (20), the module is integrated into or attached to an 
object or a device (30), electronic interface (12), and MRI system (40) with its main coil (42), three 
gradient coils (43), computer (44), gradient coils control unit (48), and image display (46). The 
MRI coils (42 and 43) are presented in more details in Figure IB with the different coils displaced 
along the MRI main axis to clarify the presentation. 

Figure 2 presents the activation sequence of the MRI gradient coils as functions of time 
during a standard spin echo scan. The steep rise and drop of the generated magnetic fields result 
with high rate of change of the magnetic flux through the coils. 

Figure 3A presents the time-derivative of the magnetic fields of the MRI gradient coils 
(which are presented in Figure 2) as functions of time. Figure 3B presents the voltages induced by 
the time-varying magnetic fields of the MRI gradient coils (those presented in Figure 2 and 3A) in 
two orthogonal sensing coils (e.g. 22, 24) as functions of time. 

Figure 4A provides a schematic configuration of three orthogonal coils (22, 24, 26) in the 
sensor (20) and the induced voltages in each coil. Figure 4B presents an example of the vectorial 
summation of the voltages induced in each coil during activation of the magnetic field of the Z- 
gradient coil into a voltage vector termed Vz . 

Figures 5 A-5D illustrate three potential configurations of coils which provide a set of three 
orthogonal coils or pairs of coils. Figure 5A shows a cubic configuration for extra-corporeal 
applications with three orthogonal coils having a typical size of up to 10 mm. Figure 5B shows a 
cubic configuration with 3 orthogonal pairs of parallel coils. Figures 5C-5D show a cylindrical 
configuration for use with catheters with a typical diameter of 2-3 mm. Figure 5C illustrates an 
axial (along the K axis) view, while Figure 5D shows a 3-dimensional display of the sensor, having 
one cylindrical coil (22) and two pairs of transverse "saddle" coils (24, 26). 
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Figure 6 presents a block diagram of the measurement and processing system, including the 
sensor (20), the electronic interface (12) and the processing and control unit (10). 

Figure 7 presents a block diagram of the preferred embodiment of the tracking methodology 
5 for various clinical uses. 

Detailed Description of the Preferred Embodiment 

Referring now to Figure 1, a typical MRJ system (40) has several modules which are 
1 0 specifically relevant to the current invention: the three gradient coils (43), the gradient coils control 
unit (48), and the image display (46). The exact implementation of the invented methodology 
depends on the MRI mode of imaging, and the following presentation relates, as a typical example, 
to a standard MRI spin-echo imaging mode. During the spin-echo protocol, repeated generation of 
magnetic fields by the 3 gradient coils provide the spatial encoding of the received MR echo and 
1 5 enable the reconstruction of the image. A sample sequence is given in Figure 2 (recorded from a 
Signa MRI system, General Electric, USA). For this sequence the system activates the Z-gradient 
coil for "slice selection", simultaneously the X and Y gradient coils for "phase encoding" and the 
X gradient coil for the "read out" phase. 

20 The gradient control unit (48) provides the processing unit (10) with real-time presentation 

of the activation sequence of the three gradient coils which generate the magnetic gradients (Figure 
2). The magnetic fields which are generated by the gradient coils have components in all three axes 
(X,Y£) y but each of the coils has a precise linear change of the amplitude of the Z-component 
along one axis only, where these coils and the generated magnetic gradients are termed by this 

25 specific axis (i.e. for the Z-gradient (Gz) the Z-component varies linearly with the Z-coordinate, 
for the X-gradient (Gx) the Z-component varies linearly with the X-coordinate, and for the Y- 
gradient (Gy) the Z-component varies linearly with the Y-coordinate). The other components of 
the magnetic fields of the gradient coils have a specific spatial distribution which depends on the 
specific design of the gradient coils. A full description of the magnetic field as function of time and 

30 location with any mode of operation (G(t,x,y,z)) can be calculated within the processing unit (10) 
by vectorial summation of the three time- variable magnetic fields of the gradient coils and the time- 
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invariant main field (Bo) of the MRI scanner (in the following presentation vectors are underlined, 
to distinguish from scalars): 

(1 ) G(t,x,y,z) = Gx(t,x,y,z) + Gy(t,x,y,z) + Gz(t,x,y,z) + Bo(x,y,z) 

5 

where x,y,z are coordinates along the three axes of the MRI coordinate system (X, Y, Z, 
respectively) and t is a time variable. Additional magnetic fields, which are generated by the RF 
(radio frequency) coils of the MRI, are not being used by the current invention. These fields, which 
alternate in the range of mega-hertz, induce high-frequency electrical potentials in the sensing coils 
1 0 which can be removed by low-pass filtration. 

In one preferred embodiment (Figure 4A), the sensor (20) consists of a set of three 
orthogonal sensing coils (22, 24, 26). The time varying magnetic field G(t,x,y,z) induces electric 
potential, or voltage (V), in each of the sensing coils, and the magnitude of the induced voltage is 
1 5 related to the time-derivative of the magnetic flux 0 through the coil, as given by Faraday Law: 

(2) V = -dG/dt 

the magnetic flux at each location is determined by the magnetic field G(t,x,y,z), the coil 
20 area (A), and the direction of the magnetic field with respect to the spatial orientation of the coil, 
as defined by a unit vector n vertical to the plane of the coil: 

(3) 0 = G(t)-nA 

25 where the dot denotes a vectorial dot product. 

Combining equations 1 - 3, the induced voltages in the coils are directly related to the time 
derivative of the magnetic field: 



30 



(4) 



V 



= -d [ (Gx(t,x,y,z) + Gy(t,x,y,z) + Gz(t,x,y,z) + Bo(x,y,z)) • n A ] /dt 
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If the sensor does not move or rotate, the Bo field and the direction vector n are constant and 
the induced voltage in each coil is given by: 

(5) V = -d [(Gx(t,x,y,z) + Gy(t,x,y,z) + Gz(t,x>y,z))] /dt • n A 

5 

The measured magnitudes of the induced voltages at the three coils and the known magnetic 
field G(t,x,y,z) as function of time at each point in the operating field (as calculated by summing 
the individual magnetic fields of all gradient coils which are active at a specific time) enable the 
estimation of the object location and direction by the following sequence of steps. This sequence 
1 0 of steps is only one option out of several possible approaches which are similar in concept and only 
differ in the actual embodiment of the estimation process. 

Step L Measurement of induced voltages 

The induced voltages in the three orthogonal coils (Figure 4) enable the calculation of the 
15 magnetic fields of the gradient coils at the location of the sensor without knowing the orientation 
of the sensor. While the magnitudes of the induced voltages at each coil change with the 
orientation, their vectorial sum is independent of the orientation and is proportional to the time- 
derivative of the magnetic field at the location of the sensor, as given by equations 4 and 5. For 
example, during activation of the Z-gradient, the time-varying magnetic field induces three voltages 
20 in the three coils. For a configuration with three orthogonal pairs of parallel coils (Figure 5B), the 
induced voltages in two parallel coils of each pair are averaged and the results are analyzed 
similarly as with three single coils. 

Thus during activation of the Z-gradient the three voltages Vzj, Vzj, Vzj^ correspond either 
25 to the measured voltages in the three single coils or are the averages of the measured voltages in 
each of the three pairs of coils. We define the induced voltages as vectors Vzj, V^j, Vz k with 
magnitudes equal to the induced voltages in each coil and directions defined by unit vectors vertical 
to the corresponding coil plane (Figure 4A). The vectorial sum of the three vectors, denoted Vz, is 
in the direction of the time-derivative of the local magnetic field of the Z-gradient: 



30 
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(6) Vz = 

-[d^tAy^ydf Oi A] a; - [d(Gz(t,x,y,z))/dfnj A] nj - [d(Gz(t,x, y> z))/dfn k A] n k 

5 This can be easily demonstrated if we break the time derivative of the magnetic field vector 

(dG(t)/dt) into three orthogonal components which are in the directions of three orthogonal coils. 
Since components parallel to the plane of each coil do not induce any voltage, the induced voltages 
VZp Vzj, Vzfc are proportional to the three components of the time derivative of the magnetic field 
and their sum is in the same direction as the time derivative of the magnetic field (dG(t)/dt). 

10 

Finally, the magnitude of the voltage vector is proportional to the magnitude of the time- 
derivative of the magnetic field of the Z-gradient at the location of the coils and at the time of 
measurement (Figure 4B): 

15 (7) |Vz| = | (d(Gz(t,x,y,z))/dt) A | 

The magnitudes and directions of the time-derivative of the local magnetic fields of the X 
and Y gradients, or of any combination of two or three magnetic fields of different gradient coils, 
are related (i.e. have the same direction and proportional magnitude) to the vectorial sum of the 
20 induced voltages in the three coils, as described above for the Z-gradient. 

The proportionality coefficient of the relation between the magnetic field and the induced 
voltage in a coil is determined by the geometry of the coils, i.e. by A, the total area of the coil (if 
a coil with multiple loops is used the total area is the sum of all areas of the individual loops). 

25 

During a typical sequence of MRI scanning two or even all the three gradient coils can be 
activated at the same time. The magnetic fields of the gradient coils are known for a specific MRI 
scanner by simulation, based on the known geometry of the gradient coils, or by measurement of 
the fields as function of location during activation of each gradient coil. The activation sequences 
30 of each gradient coil as function of time are provided by the MRI scanner as analog signals (Figure 
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2) or digital data. The known magnetic field and the activation sequence of a specific gradient coil 
can be used to calculate the magnetic field at each spatial location and for a specific time point, or 
to calculate the time derivative of the magnetic field by analog or numerical differentiation (Figure 
3A). This information can also be used to separate magnetic fields which are generated by 
5 simultaneously activated two or three gradient coils. For example, in Figure 2 the Z-gradient coil 
is activated alone, while the X-gradient coil is activated alone or together with the Y-gradient coil 
The magnitude and orientation of the magnetic field of the X-gradient coil can be determined from 
its independent activation, and this information can be used to eliminate the contribution of the 
magnetic field of the X-gradient coil from the induced voltages measured during simultaneous 
10 activation of the X and Y gradient coils and to extract the magnitude and orientation of the 
magnetic field of the Y-gradient coil. 

An alternative, more general approach is to reconstruct the reference magnetic fields, which 
are used in the estimation process (as detailed below), as a superposition of the simultaneously 

1 5 activated magnetic fields of different gradients. Thus for each time point, the activation sequences 
of the coils are used to determine the active fields and their magnitude at that time, and the overall 
field is calculated by adding the field contributions from all active coils, as shown is Equations 4 
and 5. The location of the device is estimated by comparing the measured voltages (during 
simultaneous activation of more than one gradient) to time derivative of the reference, composite 

20 magnetic field. 

Step 2: Transformation from measured voltages to magnetic fields 

The measured voltages are proportional to the time-derivative of the magnetic fields, and 

the proportionality coefficient is determined by the properties of the sensing coils (i.e. the area of 
25 each loop and the number of loops). As explained above, the time-derivative of the magnetic field 

is at the same direction as of the voltage vector (e.g. Vz for Z-gradient) and its magnitude can be 

calculated by re-arranging equation (7): 



(8) |d(Gz(t,x,y^))/dt| = |yz|/A 

30 
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Modem MRI scanners use crushers in association with each activation of a gradient. 
Typically, the crushers are spike-like rapid activation and deactivation of the gradient coil. For 
example, in a General Electric Signa MRI scanner these crushers follow the shape of a triangle 
(Figure 2) or trapezoid, and their time-derivative is similar to a positive pulse function (the up-slope 
5 of the crusher) and a negative pulse function (the down-slope of the crusher) (Figure 3A). TTie 
induced voltages are linearly related to the time derivative of the gradient field (Equations 4 and 
5) and follow the same pattern (Figure 3B). For linear activation and de-activation of the gradients, 
the induced voltages during each constant phase (i.e. up-slope and down-slope) can be averaged to 
yield a value which is directly used to calculate the amplitude of the time-derivative of the magnetic 
10 field by equation 8. Furthermore, by measuring the time of activation or de-activation of the 
gradients (e.g. At), the amplitude of the actual magnetic field can be calculated by (for linear 
activation and deactivation patterns): 

(9) Gz(t,x,y,z) = d(Gz(t,x,y,z))/dt * At 

15 

In the following presentation the determination of the location and orientation is based on 
using the magnetic field rather than their respective time-derivative. This is possible if the slope of 
the gradient activation pattern is linear and known, yet a similar procedure can be implemented by 
using the time derivatives of the magnetic fields. The magnetic fields are provided by a set of 3- 
20 dimensional maps, for example by using Cartesian coordinate system with X, Y,Z coordinates. For 
each location, the magnetic field vector can be mapped as a set of magnitude and direction 
descriptors (e.g. two angles in a spherical coordinate system), or as a set of three orthogonal 
components of the field vector. 

25 Step 3: Estimation of the location x 9 y,z of the device in the MRI coordinate system 

By knowing the 3-dimensional distributions of the magnetic fields of the X,Y and Z 
gradients (or a combination of 2 or 3 gradient fields), the instantaneous location of the device can 
be estimated. A search algorithm finds a specific location which, during the activation of the 
gradients, has magnetic fields with similar magnitudes as those calculated from the measured coil 

30 voltages. A typical search algorithm minimizes a cost function which is based on the level of 
similarity between the estimated fields and the reference known fields at the assumed location, for 
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example a least squares cost function is the sum of the squares of the differences between each of 
the estimated magnetic fields and the corresponding reference fields at the current estimated 
location. 

5 Several problems can hamper the accuracy of the estimation - the search algorithm may find 

a local minima of the cost function (i.e. a wrong solution), the cost function may be flat or noisy 
at the region of the minima which may result in a non-accurate solution, and the minimized function 
may have more than one solution (non-unique solution). 

1 0 The problem of local minima can be solved by using search algorithms which guarantee the 

convergence to the true, global minima. For example, a grid search evaluates the cost function all 
over the potential range of solutions. For the current invention, a grid-search which evaluates the 
cost function at all combinations of x,y,z coordinates at a resolution of 1 cm was found to guarantee 
convergence to the global minima. 

15 

The accuracy of the estimation critically depends on the signal-to-noise ratio of the 
measurements. When only few measurements are used, for example in our case three unknown 
location variables are calculated from only three measurements (the amplitudes of the three voltage 
vectors), any noise will bias the estimation results. The effect of noise can be reduced when more 
20 measurements are used and a least-squares estimation algorithm is applied. This can be achieved 
by using more coils, for example a set of six coils, arranged as three orthogonal parallel pairs with 
known distances between the parallel coils. Obviously, more coils will generate more data with a 
high cost of more complex processing apparatus. 

25 The problem of non-uniqueness of the solution is associated with multiple minima, for 

example due to symmetry in the cost function. The typical spatial distribution of the gradient fields 
in commercial MRI systems has symmetry in the three axes, and as a result up to eight equivalent 
minima may exist on the cost function with up to eight different solutions for the estimation 
process. Multiple solutions are a major limitation for any tracking method, and additional data must 

30 be used to reduce the number of solutions. 
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Step4: Calculation of angles between the voltage vectors 

The magnetic fields are vectors, and at each point of the imaging field the orientations of 
the magnetic fields of the three gradients are typically different, and can be used as additional 
information for the estimation process. Since the orientation of the device relatively to the MR1 
5 coordinate system is still unknown at this stage of the estimation process, the angles between the 
three gradient vectors are used instead of the global orientations of the vectors with respect to the 
MRI scanner coordinate system. The angle between any two vectors can be determined by vector 
algebra and analytic geometry* For example, the angle a between the voltage vector Vz, which is 
induced by the Z-gradient, and the voltage vector Vx, which is induced by the X-gradient, is 
10 determined by calculating the squared amplitude of the vectorial difference between the two 
vectors: 



where Vzj, Vzj Vz k and Vx^ Vxj Vx k are the measured voltages in the i j,k coils during 
activation of the Z-gradient coil and the X-gradient coil, respectively, and then calculating the angle 
between the two vectors by applying the cosine law: 



where |Vz| and |Vx| are the magnitudes of the voltage vectors induced by the Z and X 
gradients, respectively. 

The measured angles are compared to reference angular field maps, which can be generated 
25 from the 3 -dimensional field maps of the three gradients using the same approach as described by 
equations 10 and 11. 



(10) |Vz - Vx| 2 = (Vzi - Vxj) 2 + (Vzj - Vxj) 2 + (Vz k - Vx k ) : 



2 



(1 1) COS((X)= ( [Vz] 2 + |Vx| 2 - |Vz - Vx| 2 ) / (2 * |Vz| * |Vx|) 



20 



30 



In the estimation process, the measured angles are compared to the reference angles at the 
estimated location in addition to the comparison of the magnetic fields amplitudes. This additional 
information improves the accuracy of the estimation process and eliminates the problem of non- 
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uniqueness due to symmetry of the magnetic fields in the XY-plane of the MRI gradient coils. 

Using the amplitudes of the voltage vectors and the angles between the vectors, there are 
still two equivalent anti-symmetric solutions which have the same cost-function. The gradient fields 
5 of the MRI scanner are anti-symmetric - for example for a set of values of X,Y and Z coordinates 
there exists a point with the opposite X,Y and Z values (i.e. having the same absolute value but 
opposite signs) which has exactly the same absolute magnitudes and angles between the gradient 
field vectors. The distinction between the two anti-symmetric solutions can be done only during 
later stages of the estimation process, as explained below. 

10 

Following the grid search, a more accurate location can be found by local search around one 
of the two locations which were found to be the global minima of the cost-function. Since the two 
solutions are anti-symmetric, the local search can be applied around one of the two solutions and 
the final result can be used to find the anti-symmetric solution. 

15 

The local search applies a standard search algorithm, for example a Levenberg-Marquardt 
search algorithm, using either the six data points (three amplitudes of the voltage vectors and three 
angles, as detailed above), or with more data when it is available by using measurements from 
configuration with more than 3 coils. 

20 

Step 5: Determination of the device orientation 

Once the spatial location of the sensor in the magnet bore is determined through steps 1-4, 
the X,Y,Z components of the magnetic field at this location during the operation of any gradient 
or gradient combination are known for a specific MRI scanner from the reference 3 -dimensional 
25 magnetic field maps of the gradient coils. Using the voltages measured in each of the 3 coils during 
the activation of the gradients, the three rotation angles which transform from the MRI reference 
coordinate system to the local, device-attached coordinate system, are determined by an iterative 
optimization procedure. Furthermore, at this phase only one of the two anti-symmetric solutions 
provides a minimum of the new cost function, and a unique solution results. 

30 



An initial value of the three rotation angles is used to transform the X, Y,Z components of 
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the magnetic fields of the three gradients into the components of the magnetic fields in the local 
(device) coordinate system I,J,K. According to Euler's Rotation Theorem, any spatial rotation can 
be described by three rotation angles, and various conventions exist for these angles. For example 
one convention (which is typically referred to as the Euler angles), is based on rotation around the 
5 Z axis by angle "<J>", followed by rotation around the new X axis by angle "8", and finally rotation 
around the new Y axis by angle "t|r'\ The three rotations can be described by a rotation matrix: 



(12) R = 



rn ri2 ri3 

T21 T22 f23 
f31 f32 T33 



10 



15 



where the rotation matrix terms are given by: 
ri i = cos(ijr)*cos(<|))-cos(0)*sin(4))*sin(i|f); 
ri2 = cos(i|r)*sin(<|))+cos(e)*cos(<|))*sin(iJr); 
ri3 =sin(i|r)*sin(6); 

T2 1 =-sin(t|;)*cos((|>)-cos(0)*sin((|))*cos(t|r); 
T22 ~sin(i|f)*sin((J))+cos(8)*cos((|))*cos(iJf); 
r23 =cos(i|r)*sin(0); 
Oi =sin(6)*sin(4>); 
T32 =-sin(8)*cos(<|>); 
1-33 =cos(6); 



Using the rotation matrix, the magnetic field vector in the reference coordinate system of 
the MRI scanner (i.e. in the X,Y,Z system, with components Gx, Gy, Gz) can be presented in 
20 another, rotated coordinate system. If a local coordinate system I,J,K is attached to the device, and 
is rotated by the three rotation angles {4>,8,i[f } in reference to the X,Y,Z system, the three Cartesian 
components of the magnetic field vector in the rotated system (Gi, Gj, Gk) are found by: 





Gi 




m 


ri2 


ri3 




Gx 


(13) 


Gj 




T21 


T22 


T23 


• 


Gy 




Gk 




T31 


T32 


T33 




Gz 
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The calculated components of the gradient magnetic field in the I,J,K local system can be 
compared with the measured components in order to determine the three unknown angles of 
rotation. These three unknowns can be solved from the three components of one gradient field, but 
5 the results may be biased due to noise in the measurements. Better results can be achieved by using 
data from more gradients. Since all the three gradient fields are activated during every MRI 
scanning, the preferred embodiment of solving for the three angles of rotation involves the use of 
nine gradient field components, including 3 components for each of the 3 gradients fields (or 18 
components if a set of 6 coils is used) and an optimization algorithm, for example the least squares 
1 0 method which is described above, to solve for the best solution. 

Unlike the situation with absolute values of the measured voltages, which results with a 
non-unique solution composed of the true solution and an anti-symmetric one, the use of the actual 
measurements in each coil during the activation of each gradient or a combination of two or three 
15 gradients provides a unique solution. The gradient field components in the two locations, which 
correspond to the two solutions, have the same absolute values but opposite directions, so the 
induced voltages in each coil have opposite signs. Although a more rigorous mathematical analysis 
can be used to prove the uniqueness of the solution at this phase of the estimation process, a 
numerical example is provided as a simple demonstration. 

20 

For a specific location (e.g. x=20.5 cm, y=10.5 cm, z=15cm) and three rotation angles (e.g. 
<f> = -40, 0 = 80, i|r = 0) the induced voltages in the three orthogonal coils during activation of a X- 
gradient, Y-gradient and Z-gradient are given in Table 1 (units are arbitrary and the simulation is 
based on maps of the gradient fields of a Signa MRI scanner). For this location, the absolute values 

25 of the voltage vectors and the angles between these vectors are calculated and given in the Table. 
Estimation of the location, using only the amplitude of the voltage vectors, results with eight 
potential solutions which all have the same absolute values of the vectors. The angles between the 
voltage vectors are different in 6 of the 8 solutions, leaving only two equivalent solutions (the input 
location and the anti-symmetric solution x= - 20.5 cm, y = -10.5 cm, z = - 1 5cm). Comparison of 

30 the components of the voltage vectors shows that they have opposing signs in the two locations, 
which enables the determination of the true location of the sensor. 



• 
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Step 6: Improving the estimation accuracy by using measurements from all coils 
Steps 1-5 described a preferred embodiment of the invention using a two-tier estimation 
process, the first one determines the location and the second one determines the orientation of the 
object or the device, when only three orthogonal coils are used or when the measurements from the 
two coils in each pair are simply averaged. However, when all the measurements are used in the 
estimation process a more accurate estimation result can be achieved. 

The estimation process aims to find the 6 unknowns which fully define the spatial location 
and orientation of the sensor. Since the exact distance between the two coils of each pair is 
accurately known, the estimation process still aims to find 6 unknowns, for example the location 
and orientation of a set of three orthogonal coils, while the location and orientation of the second 
set of the three orthogonal coils can be defined with respect to the location and orientation of the 
first set. Thus, although we get more measurements (18 voltages for the 6 coils during operation 
of each of the three MRI gradient coils) we still have the same number of unknowns. A larger 
amount of data for the estimation process is a key for more accurate solution of the optimization 
process. 

The effect of voltages induced by the Bo field when the sensor moves or rotates 

Equation 4 provides the general description for the induced voltages in the sensor coils, but 
the description above assumes no effect of the Bo field. This assumption is correct as long as the 
sensor does not move, or when the movements are relatively slow. Since the typical rise time of 
gradients in modern MRI system is 1 millisecond (Figure 2), and body or device movements are 
typically slower (in scale of seconds or tenth of a second), the effect of the Bo can be eliminated 
by appropriate high pass filtering of the signals from the sensor coils (for example a 1 00 Hz cut-off 
frequency), and the description above can be applied on the filtered signals during movements and 
rotations of the body organs (e.g. the head) or the device. 

Yet the voltages induced in the sensor coils by the Bo magnetic field can be advantageously 
used to improve the tracking of the location and the orientation of the device or the object. Unlike 
the gradient fields which change in time, the Bo is constant and it induces electric potential in the 
sensor coils only when there is rotation of the coils which change the magnetic flux through the 
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coils. Unlike equation 3 above, the time varying variable now is the direction of the coils which is 
given by a time variable unit vector n(t): 

(14) B = G • n(t) A 

5 

By applying low-pass filter on the sensor signals the Bo-induced voltages can be extracted 
and used to estimate the time change of the sensor orientation, i.e. the three angular velocities of 
the device or the sensor. This information can be used to improve the estimation process based on 
the magnetic fields of the gradients (e.g. by providing a better initial guess to the iterative 
1 0 estimation processes) or to enable better prediction of future location and orientation of the device 
or the sensor. 

Preferred Embodiment of the Sensor 

15 The preferred, minimal configuration of the sensor includes three coils. The disclosed 

invention covers also a lesser configuration of two coils, using the same approach as presented 
above, and using the 6 potentials induced by the three MRI gradients in the two coils to calculate 
the 6 unknowns (3 locations and 3 rotation angles). However, for optimal performance more data 
should be used to reduce the effects of noise and to improve the accuracy of the tracking. A 

20 potential configuration with three orthogonal coils is presented in Figure 5 A. This configuration 
is suitable for extra-corporeal applications, for example devices for minimal invasive procedures 
like biopsy guns or surgery instruments. Furthermore, the inner space of the sensor can be used to 
contain electronic circuitry, powered by a miniature battery, for signal conditioning (e.g. filtration 
and amplification), signal transformation (e.g. into optical signal, or into frequency modulated (FM) 

25 signal), or for wireless transmission of the measured potentials. 

A more complex configuration is presented in Figure 5B, where three pairs of parallel coils 
can be used instead of three single coils, i.e., a total of 6 coils is used in a sensor. The major 
advantage of this configuration is a substantial increase in the accuracy of the tracking, since for 
30 each activation of any MRI gradient, six, rather than three, different potentials are induced, and a 
total of 18 measurements is available to estimate the 3 unknown location variables and the 3 
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unknown rotation angles with each cycle of scanning. Although the distance between each of the 
two parallel coils is small (e.g. 5-10 mm in the cubic configuration of Figure 5 A and 1-2 mm in the 
cylindrical configuration of Figure 5C and 5D), the steep gradients used with modern MRI scanners 
on one hand, and the availability of the exact distance between the two parallel coils on the other 
5 hand, enable the use of this information to increase the accuracy of the tracking. 

A second preferred configuration is presented in Figure 5C and D, and includes a cylindrical 
coil and two pairs of "saddle" coils positioned in orthogonal directions to the cylindrical coils and 
to each other (Figure 5C presents an axial view of the set of coils and Figure 5D presents an 

10 isometric view of the two pairs of saddle coils and an inner cylindrical coils, all three coils are 
axially displaced to clarify the presentation). This configuration is specifically useful for catheters 
tracking since it has a hollow cylindrical structure and it can be fixed on the tip of any catheter 
without blocking the lumen of the catheter. It can be used with stent placement apparatus, with 
various diagnostic catheters (e.g. for intra-cardiac electrophysiology studies) and with current or 

15 future therapeutic catheters (e.g. RF ablation, laser ablation, percutaneous transmyocardial 
revascularization (PMR), targeted drug delivery, local genetic substance placement, etc.). 

In a variant of the cylindrical hollow configuration the two pairs of "saddle" coils are 
replaced by two planar coils, which may be positioned inside or outside the lumen of the cylindrical 
20 coil. Although this configuration partially blocks the catheter lumen, it is simpler to manufacture 
and may be useful with applications which do not require free lumen. 

The sensors can be assembled from individual coils, for example by glueing 6 small flat 
coils on the 6 surfaces of a cube. On a catheter, one pair of coils has a cylindrical shape and can be 
25 directly wired over the shaft of the catheter, while the two other pairs have saddle shape, and can 
be glued around the cylindrical coils. Another potential approach for the construction of the multi- 
coil sensor is by using flexible printed electrical circuits, which include all the coils and are folded 
to achieve the 3-dimensional shape. 



30 
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Preferred Embodiment of the Tracking Apparatus 

The tracking apparatus (Figure 6) includes the sensor 20, the electronic interface 12, the 
processing module 1 0, and the interface with the MRI scanner. It can be custom-designed and built 
for the specific tracking application or assembled from commercially available components. 



The electronic interface (12) contains a set of amplifiers (122) to amplify the low-voltage 
potentials which are induced in the coils (from millivolts level to volts), a set of low-pass filters 
(124) to eliminate the high frequency voltages which are induced by the RF transmission, having 
frequency range of 1 0-400 MHz (depending on the MRI magnet strength), and stop-band or notch 
1 0 filter (126) to remove potentials induced by the step-wise increase of the MRI gradients, which in 
a General Electric MRI scanner produce a 128KHz artifact. Various commercial systems with 
programable amplifier\filter combinations can be used to amplify and filter the low-voltage signals 
from the sensors (e.g. SCS-802, Alligator Technologies, Costa-Mesa, CA). 

1 5 The processing and control unit (1 0) can be developed using readily available commercial 

hardware. For example, the measured signals from the sensor can be digitized by analog-to-digital 
(A/D) converter (102) using a standard data acquisition board (e.g. National Instruments, Austin, 
TX), and processed in real-time by a modem, high performance processor 104 (e.g. a Pentium III 
processor with MMX built-in DSP). Another potential solution, which provides faster estimation 

20 rates, can be based on digital signal processor (DSP) boards, having built-in or attached A/D 
converter having at least 6 channels (3 coil signals and 3 MRI gradient signals), high-performance 
DSP for iterative solution of location and orientation, sufficient memory capacity for the program 
and for data (e.g. the reference magnetic fields), and communication bus for interface with the host 
computer or directly with the MRI scanner (e.g. Blacktip-CPCI processing board and BITSI-DAQ 

25 analog input/output adapter, Bittware Research Systems, Concord, NH). The software for the DSP 
or the PC processor can be developed with standard programming languages, for example C++ or 
assembly. We have used the Matlab software development environment (The Math Works, Natick, 
MA) to rapidly implement the estimation process as described above. 



5 



30 



The interface with the MRI includes two main components - a channel to transfer the real- 
time location and orientation of the sensor, and a channel (or channels) to transfer the activation 
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pattern of the gradient coils from the MRI scanner to the processing module. Either digital 
communication channel, analog channels, or a combination of the two can be used. With the Signa 
MRI system, for example, the gradient activation sequence is available as standard analog output 
from the gradient control system, and tracking information can be received by the MRI through a 
5 standard serial communication line. 

The overall operation of the tracking system is presented below and in Figure 7. The 
induced potentials in the sensors (700), typically having a magnitude of millivolts, are amplified 
and filtered by the electronic interface module (710). The activation pattern from the MRI scanner 

10 (702) is transferred to the tracking system through the MRI interface module (704) and may be 
processed by the electronic interface module (e.g. filtered) before it is digitized by the processing 
module. The activation pattern of the MRI gradients (Figure 2) is analyzed by the processor to 
determine the activation of each of the gradient coils, e.g. by threshold triggering (712). Typically 
we will use the crushers which have longer activation times and higher amplitude of magnetic 

1 5 fields. Once an activation of any gradient coil is detected, the processor digitizes the signal from 
the coils and process it to determine the level of the induced signals (714). If the activation of the 
gradient is linear, its time-derivative during the activation is flat (Figure 3A) and the induced 
potential in the sensor coils is also flat (Figure 3B). Thus, the measured signals can be averaged as 
long as the gradient activation is on. It should be noted, however, that non-linear activation patterns 

20 can be used as long as a description of the activation pattern of the gradient coils is available. The 
measured signals from the three orthogonal coils are calibrated into magnetic fields units using the 
calibration factors of the coils (Equations 8 and 9). The measured induced voltages in the set of 
orthogonal coils are used to calculate the amplitude of the voltage-vector (Equations 6-7) and the 
angles between the voltage-vectors of the different MRI gradients (Equations 10-1 1) (716). These 

25 amplitudes and angles are used to estimate the location of the sensor in the MRI coordinate system 
(block 718, 3) and to estimate the orientation of the sensor (block 720, 4). The estimated location 
and orientation may be further processed to improve the quality of the tracking, e.g. by the 
application of a low-pass digital filter on the estimations at a specific time, using previous 
estimations, and may be transformed into a data format which is required by the MRI scanner (722). 

30 Finally, the tracking data (724) is transferred to the MRI scanner through the MRI interface module 
(704). 
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Clinical Applications 

The determined location and orientation of the sensor can be transferred to the MRI scanner 
in real-time and used for various tasks, for example for real-time control of the scanning plane, to 
display the location and orientation of the object or the device with the tracking sensor on the MR 
5 image, to correct motion artifacts. Potential clinical applications of the invention can be divided into 
applications for diagnostic MR imaging and for interventional MRI. 

Diagnostic MRI: A major problem with MR imaging is motion artifacts due to patient movement. 
With high-resolution scanning, which may require image acquisition during many seconds and even 

10 minutes, patient movement and breathing may induce motion artifacts and blurred images. MR 
scanning is specifically sensitive to movements during phase contrast angiography, diffusion 
imaging, and functional MRI with echo-planar imaging (EPI). Using the present invention for real- 
time determination of the location and orientation of the scanned object can reduce the effect of 
motion on the MR scans by real-time control and correction of the scanning plane, in order to 

1 5 compensate for the movement, or by post-acquisition image processing. 

Interventi onal MRI: The sensor can be used with various devices, like miniature tools for minimally 
invasive surgery, catheters inside blood vessels, rigid and flexible endoscopes, biopsy and aspiration 
needles. It can be used to measure the location of the device with respect to the MRI coordinate 

20 system and to enable the MR scanner to present the device location on the MR images, as visual 
feedback to the operator, or to calculate and display the line of current orientation to assist the 
operator to steer the device into a specific target. Another potential application is to slave the MRI 
plane of imaging to the tracking sensor, for example to apply high resolution imaging on a small 
volume around the site of a catheter, for better imaging of the region of interest to improve 

25 diagnostic performance or to control the effect of an intervention (e.g radio-frequency, cryo, or 
chemical ablation and laser photocoagulation can be monitored by temperature-sensitive MR 
imaging). Another potential application is to use the information of the location and orientation of 
the device in order to enable display of the MRI images in reference to the device local coordinate 
system, as if the operator is looking through the device and in the direction of the tip, similar to the 

30 use of optica] endoscopes. One more application is to use the location tracking in order to mark 
location of previous interventions on the MRI image. 
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An application with great clinical importance, where using MRI guidance is of specific 
advantage, is percutaneous myocardial revascularization (PMR). PMR is typically performed during 
cardiac catheterization. A laser transmitting catheter is inserted through the femoral artery up 
through the aorta into the left ventricle of the heart. Based on prior perfusion studies (e.g. Thallium 
5 scan) or indirect information on viability of the myocardium (e.g. by measurement of local wall 
motion), the cardiologist applies laser energy to drill miniature channels in the inner portion of the 
heart muscle, which stimulates angiogenesis and new blood vessel growth. PMR potentially 
provides a less invasive solution (compared to bypass surgery) for ischemic heart disease patients 
which cannot be adequately managed by angioplasty or stent placement. It may also be used in 

1 0 conjunction with angioplasty or stents to treat areas of the heart not completely revascularized by 
a balloon or stent placement Currently, PMR is exclusively done with X-ray guidance. The main 
advantage of MRI is the excellent performance of MRI in the assessment of myocardial blood 
perfusion, through the use of contrast agents. Thus rather than using indirect information on the 
location of poorly perfused regions, a diagnostic session of myocardial perfusion in the MRI 

1 5 scanner can be followed by immediate intervention, using the existing perfusion images and real- 
time tracking of the laser catheter with the disclosed tracking methodology. Additional advantage, 
unique to MRI, is the potential to control the intervention by high-resolution, real-time imaging of 
the myocardium during the application of the laser treatment. Furthermore, since PMR is typically 
performed on multiple locations, and a good coverage of the treated myocardium should be 

20 achieved, marking the location of the treated locations on the perfusion image, using the location 
data of the tracking system, may provide optimal coverage of the diseased region. 

Anatomically, the tracking sensor can be used for various diagnostic and interventional 
procedures inside the brain (internally through blood vessels or through burr holes in the scull), the 
25 cardiovascular system (heart chambers, coronary arteries, blood vessels), the gastro-intestinal tract 
(stomach, duodenum, biliary tract, gall bladder, intestine, colon) and the liver, the urinary system 
(bladder, ureters, kidneys), the pulmonary system (the bronchial tree or blood vessels), the skeletal 
system (joints), the reproductive tract, and others. 
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TABLE 1 

Induced voltages in 3 orthogonal coils were simulated for a sample location (X=20.5 cm, Y=10.5 cm, 
Z=l 5.0 cm) inside an MRI scanner during the activation of three gradients and are presented in the lower 
section of the table as voltage vectors (Vx, Vy, Vz). Using only the absolute amplitudes of the voltage 
5 vectors (|Vx|, |Vy|, |Vz|) results with eight different solutions (Xest, Yest, Zest). The use of the angles 
between the voltage vectors (XZjuig, YZ_ang, YXang) eliminates 6 of the solutions and provides two 
anti-symmetric equivalent solutions (1 and 8). Finally, when the three components of each of the three 
voltage vectors are used (Vx, Vy, Vz), a unique correct solution (solution 8) is obtained. 
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WHAT IS CLAIMED IS : 

1 . A method of determining the instantaneous location and orientation of an object moving 
through a three-dimensional space, comprising: 

applying to said object a coil assembly including a plurality of sensor coils having axes of 
5 known orientation with respect to each other and including components in the three orthogonal 
planes; 

generating a time-varying, three-dimensional magnetic field gradient having known 
instantaneous values of magnitude and direction; 

applying said magnetic field gradient to said space and the object moving therethrough to 
10 induce electrical potentials in said sensor coils; 

measuring the instantaneous values of said induced electrical potentials generated in said 
sensor coils; 

and processing said measured instantaneous values generated in said sensor coils, together 
with the known magnitude and direction of the generated magnetic field gradient and the known 
1 5 relative orientation of the sensor coils in said coil assembly, to compute the instantaneous location 
and orientation of said object within said space. 

2. The method according to Claim 1, wherein said magnetic field gradient is generated by 
activating gradient coils of an MRI (magnetic Resonance Imaging) scanner. 

20 

3. The method according to Claim 1 , wherein said coil assembly includes at least three sensor 
.coils oriented orthogonally with respect to each other. 

4. The method according to Claim 1 , wherein said coil assembly includes three pairs of sensor 
25 coils, in which one sensor coil in each pair has the same orientation as the other sensor coil in the 

respective pair, and in which each pair of sensor coils has a different orientation from the other pairs 
of sensor coils. 

5. The method according to Claim 4, wherein each sensor coil in a pair is parallel to, but 
30 laterally spaced from, the other sensor coil of the pair. 
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6. The method according to Claim 1 , wherein said coil assembly includes a cylindrical sensor 
coil and two pairs of sensor coils positioned orthogonally with respect to the cylindrical sensor coil. 

7. The method according to Claim 6, wherein said two pairs of sensor coils are curved and in 
5 a saddle relation to said cylindrical sensor coil. 

8. The method according to Claim 6, wherein said two pairs of sensor coils are planar. 

9. The method according to Claim 1, wherein said magnetic field gradient is generated by 
10 activating three orthogonally disposed pairs of gradient coils according to a predetermined 

activating pattern; and wherein the measured instantaneous values of said induced electrical 
potentials generated in said sensor coils are processed, together with said predetermined activating 
pattern of the gradient coils and the known relative orientation of the sensor coils, to provide an 
estimate of the location and orientation of said object. 

15 

10. The method according to Claim 1, wherein said processing is effected by an iterative 
optimization process. 



1 1 . The method according to Claim 10, wherein said iterative optimization process is effected 
20 in real time to determine the instantaneous location and orientation of said object in real time. 



12. The method according to Claim 1, wherein said method further comprises generating a 
constant magnetic field of known magnitude and direction; applying said constant magnetic field 
to said space and the object moving therethrough; and utilizing the induced electrical potentials 

25 generated in said sensor coils, upon movement of said object and sensor coils, to improve the 
computation of the instantaneous location and orientation of said object within said space. 

13. The method according to Claim 1 , wherein said object is a medical instrument moving in 
the body of a person for medical diagnostic or treatment purposes. 

30 

14. The method according to Claim 13, wherein said medical instrument is selected from the 
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group consisting of: a catheter for arteriogram, a catheter for venogram, a catheter for angioplasty, 
a catheter for stent placement, a catheter for percutaneous transmyocardial revascularization, a 
catheter for cardiac electrophysiology studies, and a catheter for gene therapy. 

5 15. The method according to Claim 1 3, wherein said medical instrument is a tool for minimal 
invasive surgery. 

16. The method according to Claim 13, wherein said medical instrument is selected from the 
group consisting of: a biopsy gun, a biopsy needle, and an aspiration needle. 

10 

17. The method according to Claim 13, wherein said medical instrument is selected from the 
group consisting of: a rigid endoscope, a flexible endoscope, a ventriculoscope, a colonoscope, a 
duodenoscope, a gastroscope, a laryngoscope, a tracheoscope, a bronchoscope, a hysteroscope, an 
urethroscope, a cystoscope, an ureteroscope, and an arthroscope. 

15 

18. The method according to Claim 1 3, wherein said medical instrument is a motion sensor. 

19. Apparatus for determining the instantaneous location and orientation of an object moving 
through a three-dimensional space, comprising: 

20 a coil assembly carried by said object and including a plurality of sensor coils having axes 

of known orientation with respect to each other and including components in the three orthogonal 
planes; 

a magnetic field generator for generating a time-varying, three-dimensional magnetic field 
gradient having known instantaneous values of magnitude and direction, in said space and the 
25 object moving therethrough to induce electrical potentials in said sensor coils; 

means for measuring the instantaneous values of said induced electrical potentials generated 
in said sensor coils; 

and a processor for processing said measured instantaneous values generated in said sensor 
coils, together with the known magnitude and direction of the generated magnetic field gradient and 
30 the known relative orientation of the sensor coils in said coil assembly, to compute the 
instantaneous location and orientation of said object within said space. 
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20. The apparatus according to Claim 19, wherein said magnetic field generator includes 
gradient coils of an MRI (Magnetic Resonance Imaging) scanner. 

21 . The apparatus according to Claim 19, wherein said coil assembly includes at least three 
5 sensor coils oriented orthogonally with respect to each other. 

22. The apparatus according to Claim 19, wherein said coil assembly includes three pairs of 
sensor coils, in which one sensor coil in each pair has the same orientation as the other sensor coil 
in the respective pair, and in which each pair of sensor coils has a different orientation from the 

1 0 other pairs of sensor coils. 

23. The apparatus according to Claim 22, wherein each sensor coil in a pair is parallel to, but 
laterally spaced from, the other sensor coil of the pair. 

15 24. The apparatus according to Claim 19, wherein said coil assembly includes a cylindrical 
sensor coil and two pairs of sensor coils positioned orthogonally with respect to the cylindrical 
sensor coil. 

25 . The apparatus according to Claim 24, wherein said two pairs of sensor coils are curved and 
20 in a saddle relation to said cylindrical sensor coil. 

26. The apparatus according to Claim 24, wherein said two pairs of sensor coils are planar. 

27. The apparatus according to Claim 1 9, wherein said magnetic field generator includes three 
25 orthogonally disposed pairs of gradient coils to be activated according to a predetermined activating 

pattern; and wherein said processor processes the measured instantaneous values of said electrical 
potentials induced in said sensor coils together with said predetermined activating pattern of the 
gradient coils and the known relative orientation of the sensor coils, to provide an estimate of the 
location and orientation of said object. 

30 

28. The apparatus according to Claim 27, wherein said processor effects an iterative 
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optimization process. 

29. The apparatus according to Claim 28, wherein said iterative optimization process is effected 
in real time to determine the instantaneous location and orientation of said object in real time. 



30. The apparatus according to Claim 1 9, wherein said magnetic field generator further includes 
a coil for generating a constant magnetic field of known magnitude and direction and for applying 
said constant magnetic field to said space and the object moving therethrough; and wherein said 
processor utilizes the electrical potentials induced in said sensor coils, upon movement of said 
1 0 sensor coils, to improve the computation of the instantaneous location and orientation of said object 
within said space. 

: 'v 31 . The apparatus according to Claim 19, wherein said object is a medical instrument moving in 
the body of a person for medical diagnostic or treatment purposes, and wherein said coil assembly 
15 is adhered to or integrated into said medical instrument. 

32 . The apparatus according to Claim 3 1 , wherein said medical instrument is selected from the 
group consisting of: a catheter for arteriogram, a catheter for venogram, a catheter for angioplasty, 
a catheter for stent placement, a catheter for percutaneous transmyocardial revascularization, a 

20 catheter for cardiac electrophysiology studies, and a catheter for gene therapy. 

33 . The apparatus according to Claim 3 1 , wherein said medical instrument is a tool for minimal 
invasive surgery. 

25 34. The apparatus according to Claim 3 1 , wherein said medical instrument is selected from the 
group consisting of: a biopsy gun, a biopsy needle, and an aspiration needle. 

35 . The apparatus according to Claim 3 1 , wherein said medical instrument is selected from the 
group consisting of: a rigid endoscope, a flexible endoscope, a ventriculoscope, a colonoscope, a 
30 duodenoscope, a gastroscope, a laryngoscope, a tracheoscope, a bronchoscope, a hysteroscope, an 
urethroscope, a cystoscope, an ureteroscope, and an arthroscope. 
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36. The apparatus according to Claim 3 1 , wherein said medical instrument is a motion sensor. 
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